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We proposed an altered configuration for dual-frequency (DP) capacitively coupled plasmas 
(CCP). In this configuration, two pairs of electrodes are arranged oppositely, and the discharg- 
ing is perpendicularly driven by two rf sources. With Particle-in-cell/Monte Carlo method, we have 
demonstrated this configuration can remove the harmful electromagnetic and DP coupling effects 
in conventional DP-CCP. 
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Dual-frequency (DF) capacitively coupled plasmas 
(CCP) are commonly used as etching and deposition de- 
vices in the microelectronics, flat panel display and solar 
cells industries [11,0]. Compare to the other two sources, 
inductive coupled plasma (ICP) and electron cyclotron 
resonance (ECR) discharges, CCP can produce uniform 
plasma over larger areas. In typical DF-CCP 0,11], the 
two rf source with different frequencies are applied to the 
same electrode or the opposite two electrodes. The high 
frequencies (hf) source controls the plasma density while 
the low frequency (If) source controls the ion flux and 
ion energy and angular distribution (lEDs and IADs). 
Quasi-independent control of plasma density and lEDFs 
is the main merit of DF-CCP over single-frequency CCP. 
Therefore DF-CCP has been widely used in state-of-art 
industry reactors. To achieve more flexibilities, some 
variants of DF-CCP, such as adopting very high fre- 
quency (VHF) 5], applying additional dc source [6} and 
series resonance CCP [7], have received intense investi- 
gation recently. 

High etching rates, uniformity, anisotropy, selectivity 
and low dielectric damage are the essential requirements 
of the CCP etching devices. There is no perfect solutions 
to satisfy all the requirements above simultaneously, one 
must make tradeoffs. Although successful in practice, 
DF-CCP still expose some problems. The first one is 
the electromagnetic (EM) effect [sl. lol. [lOj. which will oc- 
cur when very high frequency (VHF) sources (typically 
> 60MHz) are used, and therefore the excitation wave- 
length A is comparable to the reactor radius. Due to 
the standing wave effect [ll| , the merit of plasma unifor- 
mities over large areas will be broken. The second one 
is hf-lf source coupling effects When only increas- 

ing the hf voltage, the density increases and the sheath 
thickness decreases, thus more ions will appear in the 
high energy end of lEDs. This is not desired because 
dielectric damage may be introduced. But when only in- 
creasing the If voltage, the sheath thickness increases and 
the bulk plasma lengths decrease and thus the density de- 
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FIG. 1: Schematic of the perpendicularly driven DP-CCP. 



creases. This is also not desired because of the reduced 
etching rate. Some configurations, like shaped electrodes 
[l3l | and asymmetrical discharges are proposes to 
mitigate above harmful effects. 

In this letter, we propose an altered configuration for 
DF-CCP. Not like the conventions disk-like cylindrical 
configuration [l|, here the discharge chamber is a three- 
dimensional (3D) flat regular hexahedron, whose cross- 
sectional schematic is shown in Fig|T] Two pairs of 
rectangular electrodes are arranged oppositely, in the 
left /right surface and the top/bottom surface, respec- 
tively. In the front/back surface (not shown in the figure) 
and in the gaps between the electrode, thick dielectric 
insulators are placed to confine the plasmas. The dis- 
charging is perpendicularly driven by two rf sources. The 
smaller pairs (left and right) electrodes (LE and RE) are 
powered by the same hf sources, but with the phase dif- 
ference of TT. This can be realized with a opposite phase 
power divider, for example, a balun [l^ applied between 
the rf source and the devices. The bottom electrode (BE) 
is powered by the If source and the top electrode (TE) 
is grounded. The wafer is placed on the BE. We will 
demonstrate that this configuration can remove above 
problems. 

In order to study such a configuration, we have adopted 
direct implicit Particle-in-cell/Monte Carlo (PIC/MC) 
algorithms JTd\ in 2D planar geometry. Although this 
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configuration is inherently 3D, due to the confinement 
by the insulators, the plasma is uniform between the 
side- wall insulators and thus 2D model is sufhcient, if 
the front /back spacing is large. The details of the algo- 
rithms are described elsewhere (l7j , therefore we will only 
present the simulation parameters here. The BE and TE 
lengths are X = 8cm, and the LE and RE lengths are 
Y = 1.5cm, therefore we have the four gaps of 0.25cm. 
The hf sources of uJhf = 60MHz and Vhf = 50V or lOOV 
are applied to LE and RE, with the phase difference of 
TT. The If sources of ujif = 2MHz and Vif = 50V or lOOV 
are applied to BE. In the gaps, the instant potentials are 
linearly interpolated. Argon gas is used with the pres- 
sure of lOmTorr and temperature of 300K. We consider 
elastic, excitation and ionization collisions for electrons 
and elastic and charge transfer collisions for Ar"*" ions, 
respectively. Square cells are used, thus X direction is 
uniformly divided to 256 cells and Y direction divided 
to 64 cells. The space and time steps are fixed to all 
simulations. Ax = 0.02/64m, At^ = At, = 0.5 x lO'^^s. 
All results are given by averaging over one If period af- 
ter reaching equilibrium of 1000 rf periods. Because the 
diodes areas are equal for X / Y separately and only the 
ideal voltage sources in the external circuit, we did not 
consider the external circuit and the self-biasing for sim- 
pleness. 

The average potential, electron and ion density for the 
case oi Vhf = lOOV and Vif = 50V are shown in FiglH 
Near each electrode, there is a sheath, therefore the den- 
sities are in an elliptical profile and the maximum value 
is in the center. The sheathes are symmetric in both 
directions, but the average sheath thicknesses are larger 
for TE and BE than that for LE and RE. This is a nat- 
ural result, since when the density is fixed, the sheath 
thickness inversely scales with the frequency [J. 

We plotted the cross-sectional profiles of rie for differ- 
ent voltages in FiglJl The plasma density is mainly de- 
termined by the hf source. When increasing the Vhf by 
two times, the electron densities also increase by a factor 
of larger than 2. While when increasing Vf, the elec- 
tron densities are nearly unchanged. The cross-sectional 
profiles are flat in the center. As we have mentioned, 
in conventional DF-CCP, when increasing the If voltage, 
the bulk plasma length will decrease and thus the density 
will decrease [l^. In this present configuration, there is 
no such a effect, since the sheath thicknesses are decou- 
pled. 

The ion flux to BE is presented in FiglH Similar to the 
density, the flux is mainly determined by the hf source. 
When increasing the Vhf by two times, the flux also in- 
crease by a factor of larger than 2. While Vf has no 
significant effects on the flux. Compare the convention 
CCP driven by 13.56MIIz source, the flux is much larger 
even the density value is close [13]. In industry etching 
and deposition devices, it is critical to make the ion flux 
to the electrode be uniform over larger areas. As can be 
seen, over large distance (about 6cm) , the ion flux is uni- 
form. In VHF CCP, due to the finite wave lengths effect 
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FIG. 2: Average (a)electron density rie, (b)ion density rii and 
(c) potential $ for the case of 14/ = 100 V and Vif ^ 50V. 
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FIG. 3: Cross-sectional profiles of electron density Ue at Y 
1cm for different voltages. 



as we have mentioned, the plasma density and ion flux 
will not be uniform over the wafer. In this configuration, 
the wafer can be placed on X plate while the hf sources 
are in Y direction, therefore the harmful EM effects are 
naturally removed. 

In industry devices, it is also highly desired that the 
ions are anisotropic. The ion angle distributions (IADs) 
for different voltages are depicted in Figl^l Whatever 
the rf voltage is, most ions has the angle of several de- 
grees. This means the ions still keep anisotropic in this 
configuration. 

The last issue in this configuration is lEDs. The lEDs 
is influenced not only by rf frequency and voltage, but 
also by the ion transit time r^. In conventional DF-CCP, 
if one increases Vhf while keeping the Vf constant, the 
plasma density will increases and the sheath thickness 
will decrease. Therefore will decreases and the profiles 
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FIG. 4: Ion flux onto the BE for different voltages. We also 
plotted the flux to LE for comparison 
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If source and the high peak are produced by the hf source. 
For we still have uj^f > ^^i, the center energy of lEDs Vi 
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FIG. 5: Ion angle distributions (IADs) for different voltages. 
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FIG. 6: (a)Ion energy distributions (lEDs) for different po- 
sitions for the case of Vhf = lOOV and Vif — 50V; (b)Ion 
energy distributions (lEDs) for different rf voltages. 



of the lEDFs will correspondingly shift to the higher en- 
ergy tails in some cases. This will result dielectric damage 
as we have mentioned. Furthermore, the etching unifor- 
mity also requires the lEDs being unchanged at different 
position on the electrodes. 

We plotted unnormalized lEDs for different positions 
for the case of Vhf = lOOV and Vi/ = 50V in Figita). 
Here the lEDs are sampled over four 1cm segments be- 
ginning at the center of the BE. It can be seen that the 
shape of lEDs is nearly unchanged over 3cm. Note here 
the slight amplitude difference is from the flux difference 
(FigUI). Only in the electrode edge, the shape is differ- 
ent. The reason is the ions always response to the average 
electric field. In the sheath near the X electrode, the ions 
are mainly accelerated by the field Ey produced by the 
If source. 

We showed the lEDs over entire BE for different volt- 
ages in Figini^b), since the lEDs are similar at different 
position. There are clearly four peaks in all the lEDs, 
and the peaks can be divided in to two pairs. Each pair 
of the peaks seems to be produced by the hf or If source 
solely, namely, the lower peaks pair are produced by the 



still obey the simple estimation of 0.4(V// + 2Vhf), where 
2 denote the two hf sources. Unlike the conventional DF- 
CCP, the lEDs do not shift to the higher energy tail. If 
Vif >> Vfif, the lEDs will be solely determined by the If 
source. 

In summary, we have proposed an altered configura- 
tion of DF-CCP, which is in 3D flat regular hexahedron 
shape and is perpendicularly driven by two rf sources. 
The plasma density and ion flux are solely determined 
by the hf source, and are uniform over larger area, with- 
out the harmful EM and DF coupling effect. At the same 
time, the lEDs arc mainly determined by the If source, 
there are no excessive high energy ions in the tail of the 
lEDs, which will avoid the dielectric damage. In practi- 
cal devices, the geometry is 3D, but the qualitative re- 
sults here will not change and one may even adopted a 
triple frequency configuration. If X length is very large, 
rf breakdown laws by Lisovskiy should be considered 
for the reactors design. 
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